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ABSTRACT 



Context. The circumstellar environment of evolved stars is generally rich in molecular gas and dust. Typically, the entire 
environment is either oxygen-rich or carbon-rich, depending on the evolution of the central star. 

Aims. In this paper we discuss three evolved disc sources with evidence of atypical emission lines in their infrared spectra. 

The stars were taken from a larger sample of post-AGB binaries for which we have Spitzer infrared spectra, characterised 

by the presence of a stable oxygen-rich circumbinary disc. Our previous studies have shown that the infrared spectra 

of post-AGB disc sources are dominated by silicate dust emission, often with an extremely high crystallinity fraction. 

However, the three sources described here are selected because they show a peculiar molecular chemistry. 

Methods. Using Spitzer infrared spectroscopy, we study in detail the peculiar mineralogy of the three sample stars. 

Using the observed emission features, we identify the different observed dust, molecular and gas species. 

Results. The infrared spectra show emission features due to various oxygen-rich dust components, as well as CO2 gas. 

All three sources show the strong infrared bands generally ascribed to polycyclic aromatic hydrocarbons. Furthermore, 

two sample sources show Ceo fuUerene bands. 

Conclusions. Even though the majority of post-AGB disc sources are dominated by silicate dust in their circumstellar 
environment, we do find evidence that, for some sources at least, additional processing must occur to explain the 
presence of large carbonaceous molecules. There is evidence that some of these sources are still oxygen-rich, which 
makes the detection of these molecules even more surprising. 
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1. Introduction 

The chemistry of the circumstellar environment (CE) of 
evolved stars can typically be described as either oxygen- 
rich or carbon-rich, and reflects the evolution of the central 
object along the Hertzsprung- Russell diagram. Stars on the 
Asymptotic Giant Branch (AGB) start as oxygen-rich stars, 
but their carbon abundance rises during the evolution on 
the AGB, since the carbon produced in the core is dredged- 
' up to the stellar photosphere. If this dredge-up process re- 
peats enough times, carbon becomes more abundant than 
oxygen: the star has become a carbon star. The main pa- 
rameters determining the mass of carbon dredged up are 
the initial m ass of the central star and its metallicity (see 
lHerwidF2005[ for a detailed description of AGB evolution). 



* This research has made use of the SIMBAD database, oper- 
ated at CDS, Strasbourg, France. Based on observations made 
with the SPITZER Space Telescope (program id 3274, 50092), 
which is operated by the Jet Propulsion Laboratory, California 
Institute of Technology under a contract with NASA. 
** Postdoctoral Fellow of the Fund for Scientific Research, 
Flanders 



Depending on the C/0 ratio of the central star, its CE 
will be either oxygen-rich or carbon-rich. The less abun- 
dant species will be locked in carbon monoxide (CO), 
which is an extremely stable molecule and one of the 
first to be formed. If the central star is oxygen-rich, the 
GE is typically characterised by gas and dust features 



such as from OH, GO2, HoO, oxides, and silicates 



'Ju sttaiiont et all Il996l: IWaters et al.l 119961: iMolster et cdl 
i2002al : iGielen et al.l 12011 ). For some oxygen- rich sources, 
the presence of carbon-rich molecules ca n be explained 
by no n-equilibrium effects, such as shocks (^ Neiad fc Millail 
1 1981 iDuari et al.l [l999: £herchncfl 2006). For example, 
S-type stars, which are undergoing the transition from 
oxygen-rich to carbon- rich and thus have C/0^ 1, do 
not yet form carbon-rich dust. However, due to non- 
equilibrium chemistry effects carbon-ric h molecules , such 
as PAHs, HCN , and C2H2 can be seen (|Honv et al.l 120091: 
ISmolders et al.l l2010t ). In the carbon-rich case, the ex- 
cess carbon is used to form, for example, GN, G2, G2H2, 
CH4, polycyclic aromati c hydrocarboiis (PA Hs), and SiG. 



(e.g. Bakkcr ct al. 1997; Speck et al. 1997; Pcctcrs etal 



12002. : .Honv et al., ,2003^ ,Matsuura et al.. ,2006,: ,Volk et al. 
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l201l|) . Recently, even more complex carbon molecules, 
namely Cgo and C70 fuUerencs, were detected in the CE 
of evolved s tars (Cami ct al. 2010; Garcfa-Her nandez et al.l 
I2OIOI 1201 11: IZhang fc Kwodi201lll . 

In some cases, a stron g mixed chemistry is seen in 
the CE of evolved objects (jSzczerba et al.l [20071 ) ■ One of 
these is a particular class of AGB stars, the silicate car- 
bon stars, where the infrared spectr a show features of both 
C-rich species and O-rich silicates (lYamamura et"aI1l2000t 
iMolster et aU I2OOII: iGarda-Hernandez et al.l I2006D . Also 
some post-AGB stars and proto-planetary nebulae (pPNe) 
show features of both PAHs and silicates in their spectra 
(iBeintema et al.lll996l: IWaters et al.lll998l: iMatsuura et a,l.l 



12004 ICerrigone et al.fl2009HGuzman-Ramirez et al.ll2011h 

For all these objects the dual chemistry is usually a re- 
sult of the evolution of the central star. The O-rich mat- 
ter was expelled when the star was still oxygen-rich, fol- 
lowed by a transition from C/0< 1 to C/0> 1 in the stel- 
lar photosphere. Subsequent outflow of material will then 
be carbon-rich. This effect is even more pronounced if the 
O-rich material is stored in a stable disc around the cen- 
tral star, as is the case for the silicat e carbon stars and 
several of the mixed-chemistrv pPNe ([Waters et al.l Il998t 
iDeroo et al.]|2007t ICerrigone et al.ll2009ll 

In this paper we look in more detail to a sample of three 
post-AGB disc sources. These sources are part of the larger 
sample of 57 post-AG B disc source s with Spi tzer infrared 
spectr a, as discussed in lGielen et al.l (|201ll ). In lGielen et al.l 
([20T1I) we find that the circumstellar environment of these 
post-AGB disc sources is dominated by emission of silicate 
dust. However, four sources do show features attributed 
to carbonaceous molecules, in the form of PAHs and Ceo 
fuUerenes, besides the silicate emission. The carbonaceous 
content of one of th ese sources has already been discussed in 
iGielen et al.l (|2009D . In this paper we focus on the remaining 
three stars. 

Emission of PAHs is alread5;_ observed in 



sev- 



eral po st-AGB sources ( e.g. Bein tema et al.l 19961: 



Cami fc Y amamura 2001: Pce ters et al.i2 002: Molster et al.l 



Matsuura et al.l 120041: 



Gielen et al.l 120091 ). where 
the central stars appear carbon-rich. There 



IS 



2002bl: 

most of 

evidence that the central stars from our sample are oxygen- 
rich, which would make the detection of carbonaceous 
molecules even more surprising. Moreover, this is the first 
time fullerenes are detected in O-rich post-AGB sources. 

The outline of the paper is as follows: in Section [2] we 
present the programme stars and observation details. All 
sample stars show features due to O-rich dust species, which 
are discussed in Sect.[3| O-rich CO2 gas emission is also seen 
in a few sources, and discussed in Sect.^j Clear emission of 
PAHs can be seen in the 6 — 12 /zm region of the observed 
infrared spectra, as can be seen in Sect.[5l In Sect.[6]we re- 
port on the detection of the Cgo fullerenes and derive some 
basic parameters such as the excitation temperature and 
mass. We end with discussions and conclusions in Sects. [7| 
andm 

2. Programme stars and observations 

The binary post-AGB stars discussed here are those sources 
from t he larger Spitzer s ample of post-AGB disc candi- 
dates ([Gielen et al.l l201l[ ). showing features of carbona- 
ceous molecules besides the dominant silicate emission 
seen in these stars. The sample of disc candidates was 



selected on the basis of their infrar ed colours, a s dis- 
cussed in iDe Ruvter et al.l ([2006[ ) and iGielen et all ([20 lit) 
These binary post-AGB sources are characterised by the 
presence of a stable O-rich dust disc, with strong fea- 
tures of am orphous and crystalline silicates ([Gielen et al.l 
I2OO8I [201 ID . These post-AGB disc sources typically have 
an O-rich stellar photosphere, with no evidence for 
strong third dredge- up on the AGB, as determined, 
for example, by the very low C/0 ratio (<< 1), low 
i2Qy'i3(-; ratio, and the lack of s-process enhancemen t (e.g . 
Gonzalez et^ll997l: [V an Win ckel 200l iMaas et al.ll2005l: 
Revniers fc Van Winckeil 120071: ICielen et al.ll2009t) . This is 



surprising since several stars have initial ma sses or lumi- 
nosit i es hig h enough to evolve to a carbon star ([Gielen et al.l 
[20091 120111 ). Apparently, the AGB evolution of these stars 
was shortcut, probably under the influence of strong binary 
interaction. 

All stars were observed with the Spitzer-IRS spectro- 
graph, in high- or low-re solution modes (resolution of ^ 600 
and ^ 160 respectively [Werner et al.l l2004at iHouck et al.l 
^2004). Exposure times were chosen to achieve a S/N of 
- 400, and resulted in values of 2200 s, 1750 s, and 300 s 
for IRAS 13258, IRAS 06338, and HD 52961, respectively. 
The spectra were extracted from the droop and RSC data 
products, processed through the S18.0 version of the SSC 
(Spitzer Science Centre) data pipeline. For the spectral ex- 
traction we used data reduction packages, developed for 
the c2d (Core to discs) and feps (Formation and evolution 
of planetary systems) key projects. For a det ailed descrip- 
tion of these reduction packages, we refer to iLahuis et al.l 
([2006D and lHines et~an ([2005[ ). The resuhing infrared spec- 
tr a for our sam p le sta rs can be seen in Figure [T] As shown 
in IGielen et al.l ([20111 ). the infrared spectra of these three 
stars show emission features which can not be explained by 
the typical silicate dust species, as is the case for the rest 
of the larger sample post-AGB disc sources. Details for the 
programme stars can be found in Table [1] 

For IRAS 06338-1-5333 (or HD 46703, IRAS 06338 here- 
after) and HD 52961 the binarity is conflrmed, and or- 
bital parameters are given in Table [T] These two stars also 
show the depletion patt ern of condesable elements in the 
stella r photosphere (Waelkens et al.l Il991bl : iHrivnak et al.l 
'2OO 8D, commonly seen in post- AGB disc sources (e.g. 
Gi ridhar et al.l [20051: iMaas et al.l [2005iV This abundance 
pattern, characterised by a lack of refractory elements, is 
the result of gas-dust separation in the disc, follow ed by a 
re-accretion of the cleaned gas component (Wa ters et al.l 
1992). HD 52961 is an extreme example of this pro- 
cess, having a [Fe/H] ratio of —4.8 and [Zn/Fe] = -1-3.1. 
Chemical studies have shown that the sta r s are oxygen- 
rich, with C/0 < 1 (iLuck fc Bond [1 981 'Bond fc Lucd 
[19871: IWaelk e ns et al.l Il991bl: IVanW inckcl et al. 1992 
b akker et ati Il997t iHrivnak et al.l 1200 8). Unfortunately, 
no C/0 ratio has been determined for IRAS 13258-8103 
(IRAS 13258 hereafter). Optical spectra obtained so far 
cannot give conclusive evidence for either a carbon- or 
oxygen-rich nature of the stellar photosphere (Thomas 
Lloyd Evans, private communication). So far, all post- 
AGB disc sources where we have stu died the stellar photo- 
sphere, show an O-rich chemistrv fe.g.lGonzalez et alll997l : 



iGiridhar et al]ll998l 120001: IVan Winckei2003l) . So it is most 
likely that this star also follows this trend. If not, and the 
photosphere is C-rich in stead, this would be the first C-rich 
post-AGB disc source detected. 
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Table 1. The name, equatorial coordinates a and (5(J2000), spectral type, effec tive temperature Teg, su r face gravity logg 
and mctallic ity [Fe/H ] of our sample stars. For the model parameters we refer to lWaelkens et all (|l991bf ): |De Ruvter et al.l 
(2006); Hri vnak et al., (2008). The spectral type for IRAS 13258 was determined from optical s pectroscopy by Thomas 
Lloyd Evans (pri vate communicatio n). Also given is the orbital period ( see references in iDe Ruv ter et al.| 20061: 



Gielen et al.]|2007D a.nd th e C/0 ratio (|Luck fc Bo^ll984l:[Bo"nd k Luck|[T987l: IWaelkens et aLllTOQlhiVan Winckel et al.l 
1992t iHrivnak et al.l[2008l) 



Name 


Q (J2000) 


5 (J2Q0Q) 






lFe/H| 




C/0 




(h m s) 




(K) 


(cgs) 




(days) 




IRAS 06338+5333 


06 37 52.4 


+53 31 02 


6250 


1.0 


-1.6 


600 


0.7 


IRAS 13258-8103 


13 31 07.1 


-81 18 30 


F4Ib-G0Ib 










HD 52961 


07 03 39.6 


+10 46 13 


6000 


0.5 


-4.8 


1310 


0.7 



~i — ' — ' — ' — ' — r 




3, 0.3 



15 20 25 

Wavelength [micron] 

Fig. 1. The Spitzer-IRS high- and low-resolution spectra 
of our three sample stars IRAS 13258, IRAS 06338 and 
HD 52961, together with the spectrum of EPLyr, which 
shows similar emission features. The spectra are scaled and 
offset for comparison. The inset plot shows the 7 fj,m region 
of IRAS 06338 and IRAS 13258. The dotted gray lines mark 
the observed features at 6.63, 6.83 and 7.07 /xm, probably 
due to small PAH, aliphatic hydrocarbons, and/or hydro- 
genated amorphous carbon emission. 




6 8 10 12 

Wavelength [micron] 



6 8 10 12 
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Fig. 2. Results of our modelling of the 10/ini feature 
in IRAS 06338, using silica (model A, left) and alumino- 
silicates (model B, right). The different dust contribu- 
tions are: amorphous pyroxene (violet), amorphous silica 
(blue), enstatite (magenta), forsterite (light green) and 
NaAlSi40io (cyan). The combined spectrum is given in red. 
On the bottom we show the residuals of this model, and 
overplot class B (orange) and cl ass C (dark green) PAH 
emission feat ures as described bv lPeeters et al.l (|2002[ ) and 
ivan Dicdcnhoven et al.l (|2004l) . 



All three sample stars show the typical broad infrared 
excess, associated with the presence of a circu mbinary disc 
(|De Ruvter et al.ll2006t iGiridhar et alll2005D . This excess 
stars already at near-infrared wavelengths, and is very dif- 
ferent from the excess related to an outflowing shell, where 
the excess starts and peaks at much longer wavelengths. 
However, HD 52961 and IRAS 06338 do not show the very 
high Lifj lLtf ratio, seen in o ther post-AGB disc sources (see 
Fig. 1 in lGielen et aI1l201iri . For HD 52961 Lm/L^ = 12% 
and IRAS 06338 has a ratio of only 3%, whe reas the typ- 
ical value for such sources is ar ound 40-50% (|Gielen et al.l 
l2008ll201lD . Ahhough for IRAS 13258 the lack of stellar pa- 
rameters and value for the total reddening makes an exact 
determination of the Lm/ L^, ratio difficult, th e SED shows 
that t he infrared ex c ess is significant (> 50%. lGielen et al.l 
120111 ). iLlovd Evani ()1997f ) deduced that for this star the 
disc is most likely seen very edge on, since the strong Na D 
emission lines in the optical spectrum show that the star is 
seen by reflection on the circumstellar matter, which points 
to a close to edge-on orientation of the disc. 



IRAS 06338 and HD 52961 have strong similarities with 
another post-AGB disc source, namely EPLyr. This stars 
also has an O-rich classification (jGonzalez et al.iri997( ) , and 
shows strong class C PAH features and CO2 gas line emis- 
sion (see Fig. HI), together with stro ng features due to crys- 
talline silicates (jGielen et al.ll20d9h . Furthermore, EPLyr 
also has an atypically low Ljp/L^, ratio of 3%. 

3. Silicate dust features 

In order to study the carbonaceous molecules visible in the 
spectra, we need to remove the underlying silicate dust con- 
tribution. Our analysis of the residual carbonaceous emis- 
sion is then further discussed in Sections |S| and jS] 

3.1. The lOfim region 

The spectrum of IRAS 06338 is dominated by a broad 
emission feature peaking at 9.2 (see Fig. [2]). The 
emission is clearly different from the typical 9.8 fim fea- 
ture seen in other post-AGB disc sources (see Fig. 6 in 
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iGielen et al.ll20li[) , where the feature is mainly due to emis- 
sion of amorphous and crystalline olivine. Clearly, other 
dust species, which peak at shorter wavelengths, are needed 
to explain the observed emission. One possibility is sil- 
ica, w hich has different p olymorphs that all peak around 
9/im (jSargent et al.|[2009l ). To determine the different con- 
tributing dust species to this feature, we assume the emis- 
sion to come from an optically thin region and linearly 
add d ust absorption profiles, as done in I Gielen et al.l ()2008l 
1201111 . The mass absorption coefficients for the dust species 
are calculated from refractory indice s in gaussian random 
fields (GRF) dust approximation fShkuratov fc Grvnkol 
|2005|) . The refra ctive indices for the different dust species 



are ta ke n from iServoin fc Piroul (Il973l) : iDorschner et al 
19951) : iHenning fc Mutschkd (|1997[ ). and iJaeger et al 



1998f) . Detailed results of the modelling can be found in 

Table 

We can reproduce the feature using amorphous pyrox- 
ene, amorphous silica and forsterite (Fig. [2] Left). The 
residual emission at 6.04 — 6.28 — 8.23 and 11.28 fiui can be 
explained by emission due to PAHs. Note that we do not 
find evidence for the presence of amorphous olivine, even 
though it is usually one of the most ab undant sp ecies seen 
in the spectra of post-AGB disc sources (iGielen e t al. 2008., 
1201111 . and we do see features due to crystalline olivine at 
longer wavelengths (see Sect. 13.21). Sonie alum ino-silicates 
also peak around 9 /zm ([Mutschke et al.lll998D . and a good 
fit can also be obtained using NaAlSi40io, forsterite and en- 
statite (Fig. O Right). The residual spectrum is then very 
similar to the former case, but the 8.23 /im residual feature 
is now blueshifted to 8.00 ^m, with an additional feature 
at 8.66 /im. Again, the residual features can be explained 
by emission due to PAHs. In none of the models, adding 
amorphous olivine improved the fit. The models described 
here are most likely not the only solutions, since emission 
features of amorphous dust are often relatively interchange- 
able. However, all models show strong residual emission at 
8 /im, but the exact shape remains unclear. 

HD 52961 also shows a prominent feature in the 10 fj,m 
region, peaking around 9.7 /im, pointing to the presence 
of amorphous olivine and pyroxene. Again we see a clear 
shoulder at 8 /.tm and a strong feature at 11.3 /im. This fea- 
ture is much broader than what is observed for IRAS 06338 
and is most likely a blend of forsterite and PAH emis- 
sion. We are able to get a reasonable fit to the spectrum 
using amorphous olivine/ pyroxene, amorphous silica and 
forsterite (see Fig. [3]). However, it proves difficult to re- 
produce the strong 11.3 ^m without getting too much flux 
in the 10.5 fim region. This could be due to the adopted 
synthetic spectrum for forsterite. The synthetic spectra of 
crystalline silicates are known to be very sensitive to the 
exact chemical co mposition, t emperature, the grain shape 
and grain size fe.g. lKoike eta l. 2003, 2010). Similarly, also 
the very strong observed 16 /im feature, generally associ- 
ated with forsterite emission, proved dif ficult to reproduce 
with current forsterite synthetic spectra (jGielen et al.ll2009l 
I2011h . As for IRAS 06338, we again find residual emission 
at wavelengths corresponding to PAH emission features. 

IRAS 13258 shows no strong silicate features in the 
10 fj-Tn region. A weak feature can be seen ranging from 
8 to 13 ^m. This could be due to weak emission of amor- 
phous silicates. Since there is evidence that we see the disc 
close to edge on, this could be due to heavy reddening of 



3, 3 




9 10 11 

Wavelength [micron] 

Fig. 3. Results of our modelling of the 10 fj,m feature in 
HD 52961. The different dust contributions are: amorphous 
pyroxene (violet), amorphous olivine (turquoise), amor- 
phous silica (blue) and forsterite (light green). The com- 
bined spectrum is given in red, with the continuum plotted 
in gray. On the bottom we show the scaled residuals of this 
model, and overplot an example of class B (orange) and 
class C (dark green) PAH emission features. 

the central star and the hot inner regions of the disc, where 
one would expect most of the 10 fim emission to arise from. 

3.2. The 18-36 fitn region 

For all three stars, the observed long-wavelength spectra 
show emission features that can be reproduced by emis- 
sion of amorphous olivine/pyroxene, forsterite/enstatite 
and amorphous silica. In Figure |4] we plot the observed 
long-wavelength spectra, together with model spectra con- 
sisting of the above mentioned dust species. Detailed re- 

Since 



suits of the modelling can be found in Table lA.ll 
the limited wavelength range makes it difficult to distin- 
guish the underlying continuum from contribution of amor- 
phous olivine/pyroxene, we only include amorphous sil- 
ica, forsterite and enstatite in our modelling. NaAlSi40io, 
which we used in the 10 ^m modeUing of IRAS 06338, has a 
similar feature as silica at 21 /xm, and so both dust species 
are interchangeable in this model. The features at 21 /im of 
NaAlSi40io and Si02 are relatively weak compared to the 
10 /im feature, especially at higher temperatures, so we do 
not expect to see a very strong contribution of these species 
at 21 fim. 

For IRAS 06338 and HD 52961, forsterite is clearly the 
dominant crystalline species, with a smaller contribution of 
enstatite and amorhous silica. In IRAS 13258, there is a sig- 
nificant contribution of enstatite, but there is no evidence 
for amorphous silica. For IRAS 13258, the strong forsterite 
33.6 /im feature points to the presence of a large fraction of 
cool crystalline silicates in the disc. If we see the disc very 
edge on, we might see more of the cooler outlying mate- 
rial, explaining the dominance of the cool long-wavelength 
features in the spectrum. 

4. CO2 

In IRAS 06338 and HD 52961 we see strong narrow lines in 
the 13 — 17 /im region, which can be identified with CO2 
gas emission (Fig. [S|). 
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20 25 30 35 

Wavelength [micron] 



Fig. 4. The long-wavelength Spitzer spectrum of our sam- 
ple stars, together with the model spectrum (red) consisting 
of amorphous silica (blue), forsterite (light green) and en- 
statite (magenta). The adopted continuum is given in gray. 
The spectra are scaled and offset for comparison. 



could contribute to the observed feature in IRAS 06338, but 
additional emission is still needed to give rise to such broad 
emission. Only a few dust species have strong emission fea- 
tures in the 15 fim region, most of them are titanium-oxide 
combination, but none of those species provides an expla- 
nation for the observed feature. 

A likely explanation for the origin of this broad emis- 
sion band is that there is a large column of optically thick 
CO2 gas. Simple isothermal LTE models of optically thick 
CO2 gas show a broad and strong plateau underneath the 
strong Q-branch bands. This CO2 plateau is made up by 
the multitudes of ro- vibrational lines that bunch up in this 
wavelength range. Such plateaus can indeed reproduce the 
overall width of the observed emission plateau, but overes- 
timate the flux shortward of 15 ^m. However, if the gas is 
optically thick, the spectrum can be very sensitive to tem- 
perature and density gradients in the gas, and it is thus 
conceivable that a more adequate model would in fact re- 
produce the observations quite well. Such is also the case 
for HR4049, where optically thick models can simultane- 
ously explain the broad emission plateau and the individual 
strong Q-branch bands (Malek et al., in preparation). 




13 14 15 16 17 18 

Wavelength [micron] 



Fig. 5. The CO2 emission region of IRAS 06338, com- 
pared to_othe£_20st;^AGBsources^^ emis- 
sion (|Cami fc Yamamural[200TI: iGielen et alll2009h . Dotted 
lines show the location of ^^C^^02 features, dashed lines 
give "CI6O2 and dot-dashed lines show ^'^O'^^C^^O. Dot- 
dot-dashed lines mark unidentified features. The spectra 
are scaled and offset for comparison. For HD 52961, the un- 
derlying feature at 16 /im was removed. 

For HD 52961, an underlying feature at 16 /im is seen, 
probably (partly) coming from crystalline silicates, more 
specifically forsterite. The CO2 gas emission in IRAS 06338 
is similar to the e mission seen in HD 52961 and EP Lyr 
iGielen et al.l()2009l ). and can also be identified with the CO2 
isotopes 13C1602 and Unfortunately, 

the low spectral resolution of IRAS 06338 does not allow to 
determine the exact isotope abundances, as was done for 
HD 52961. 

The CO2 lines in IRAS 06388 seem to by lying on top 
a very broad (~ 4 /im ) feature at 15 /im. Note that the 
15 /im feature observed in IRAS 06338 is much broader than 
the 16/im feature we see in HD 52961, and thus probably 
does not have the same origin. Of course, this 16 /im feature 



5. PAHs 

After the removal of the contribution of the O-rich dust 
(Sect. [3]), we find clear evidence for PAH features in all 
sample stars. The PAH emission features are generally di- 
vided in three classes, dependin g on the central wav elength 
of the most prominent features (jPeeters et al.|[200^ . Class 
A sources have features at 6.22, 7.6 and 8.6 /im. Class B 
sources show quite some variability in peak positions, but 
tend to be shifted more to the red than class A features. 
The more rare class C sources show emission features at 
6.3 /im, no emission near 7.6 /im, and a broad feature cen- 
tred around 8.2 /im, extending beyond 9 /im. 

As already discussed in Section [31 the exact identifica- 
tion of the PAH emission in the 8 /im region for IRAS 06338 
is uncertain. The features at 6.04 — 6.28 — 11.28 and 12.6 /tm 
are relatively independent from the underlying dust emis- 
sion, unlike the 8 /tm feature, where the peak position can 
shift from 8 to 8.3 /im, depending on the composition of 
the underlying silicate dust (see Fig. [2]). When peaking at 
8/im, the feature resembles more a class B PAH feature. 
With the peak at 8.3 /im, the feature is identified as class 
C. 

IRAS 06338 also shows very unusual emission in the 
7/im region. Three clear features at 6.63, 6.83 and 7.07 /im 
can be (see inset Fig. [T]), together with two smaller fea- 
tures at 6.49 and 7.25 /im. The features are re miniscent of 
the nar row features seen in the PAH spectra of ICami et al.l 
(j2011bL Fig. 2), where the features are due small PAHs, 
consisting of less than 30 carbon atoms. Similar features at 
6.85 and 7.25 /im have been observed in a few other sources 
and identified with aliphati c hydrocarbons and/or hydro- 
genatcd amorphous carbon (jFurton et al.lll999l : IChiar et al.l 
2000; Sloan et al. 200^. A smah feature at 6. 6-6. 7/im has 
been observed in other circumstellar en vironments, and ten- 
tatively identified w i th PAH emission (jPeeters et al.l 119991 : 
iWerner et af]l2004bl: IS^th et al.l[200l . 

The spectrum of IRAS 13258, shortward of 15 /tm, is 
dominated by PAH emission, with features at 6.02, 6.30, 
8.22 and 11.30 /tm. The large PAH feature around 8/im is 
again relatively broad, and appears to be a blend of smaller 
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Fig. 6. Identification of tlie PAH features in IRAS 13258. 
On the top we plot the spectrum shortwards of ISfiui, to- 
gether with a spline fit to represent the underlying con- 
tinuum. The residual spectrum is given below, together 
with an example of the class B (orang e) and class C (dark 
gree n) PAH features as described by f eeters et al.l (|2002[ ) 
and Ivan Diedenhoven et al.l (|2004ll . 



features, with peaks at 7.65, 8.22 and 8.58 /im. This broad 
band resembles the class C profile but is clearly narrower 
and has a red wing that resembles more a highly redshifted 
B profile (peaking at 8 //m) . Given that the peak position is 
8.22 fim, we classify it as a C profile. Other small peaks can 
be seen at 9.1, 9.6, 10.0, 10.6 and 12.6/im. The spectrum 
shows strong similarities with that of EP Lyr (see Fig. [1]) , 
which also shows PAHs emission at shorter wavelengths and 
strong crystahine silicate emission at longer wavelengths 
(jGielen et al.l[2009l) . 

For HD 52961 the residual spectrum (Fig. [3]) shows sim- 
ilar PAH features as for the previous two sources. We again 
find a broad 8 /im feature, with a possible small side feature 
at 9.2 /im . This would then make it a BC class (Sloa n et al.l 
12005 1 The 11.3 /im feature is very similar to that observed 
in IRAS 13258, with a similar width and central wavelength. 



6. Fullerenes: Ceo 

As can be seen already from Figure. [U the spectra of 
HD 52961 and IRAS 06338 show weak features near 17.4 
and 19 /im (Fig [7]). These features have recently been iden- 
tified with the infrared active vibrational modes of neutral 
Ceo in the Spitzer spectr um of the C-rich planetary nebula 
Tc 1 (jCami et al.ll2010l l. In that source, the bands show 
fairly broad emission bands with a roughly Gaussian pro- 
file, similar as in our objects. However, in our spectra, the 
short-wavelength bands at 7 and 8.5 /im are missing, as are 
bands of C70. 

The two Cgo features are weak, and especially the 
17.4 /im band is furthermore contaminated by some CO2 
emission (see Sect. H]). We can get an estimate of the phys- 
ical conditions and ma sses involved , follow ing a similar ap- 
proach as described in ICami et al] ()2010() , using the emit- 
ted power in each of the bands. Given the absence of the 
bands at shorter wavelengths and contamination with CO2 
lines (see Sect. [5]), the diagnostic value of the Cgo bands is 
however somewhat limited. 



6.1. HD 52961 

The band profiles of the Cgo bands in HD 52961 appear 
fairly symmetric, and Figure [7] shows that a gaussian pro- 
file reproduces the observations well. From integrating the 
spectrum, we find that the total power emitted in the Cgo 
bands is 2.6±0.2 x lO'^^ and 6.5±0.2 x lO"^*' Wm'^ for the 
17.4 and 18.9 /im bands respectively. From the ratio of these 
numbers, we then find a nominal excitation temperature of 
152^22 K . At a distance of 2.1 kpc, the observed emission 
requires a total mass of about 3.3 x 10^* M©. Note that 
these uncertainties do not include estimates of systematic 
errors made by propositioning of the continuum level. 

The observed widths of the features in HD 52961 are 
some what narrower tha n in the C-rich planetary nebula 
Tc 1 (jCami et al.ll2010l ). Smaller widths are expected for 
lower temperatures, and indeed, also the excitation tem- 
perature for HD 52961 is somewhat lower than the ^330K 
found for Tc 1. Note though that the observed widths are 
smaller than what would be expected based on labora- 
tory experiments where band widths have been measured 
at various temperatures (e.g. iNemes et al.l 11994 ). The low 
temperature also explains the absence of the two bands at 
shorter wavelengths: at an excitation temperature of 200 K 
and assuming thermal excitation, the total power in each 
of the bands at shorter wavelengths is less than 10% of the 
power in the 17.4 /im band. 




16.5 17.0 17.5 18.0 18.5 19.0 19.5 20.0 
Wavelength [/im] 



Fig. 7. The observed spectra of HD 52961 and IRAS 06338 
(black) with the estimated spline continuum (grey). In red 
we overplot the Gaussian profiles that were fit to the two 
Cgo bands. The spectra are scaled and offset for compari- 
son. 



6.2. IRAS 06338 

The IRAS 06338 observations correspond to a much lower 
spectral resolution than those of HD 52961, and at the 
same time the Cgo features are weaker, and heavily con- 
taminated by high-excitation emission bands due to CO2. 
This makes a good determination of the relevant param- 
eters much harder, and any derived quantity is probably 
rather uncertain given the systematic errors possible on the 
local continuum. We masked the CO2 bands, and adopted 
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a local continuum, and fitted Gaussian profiles to the bands 
(See Fig. [21). 

The widths (FWHM) of the 17.4 and 18.9 ^im bands 
are comparable to what we find for HD 52961. The emitted 
power is 2.7 x 10"^'' and 7.0 x 10~^''Wm~^ respectively, 
yielding a nominal excitation temperature of 145 K and a 
Cgo mass of 1.5 x 10^* (adopting a distance of 3.9 kpc). 

7. Discussion 

7.1. Mineralogy 

In all three sample sources, we find evidence for O-rich sil- 
icate dust features. Both IRAS 06338 and HD 52961 have a 
peculiar silicate mineralogy, very diffe rent to what is gen - 
erally seen in post-AGB disc sources llGielen et al.ll201lD . 
The sources seem to be dominated by warmer dust species, 
with strong emission in the 10 /im region, but very weak 
features at longer wavelengths. This is also seen in interfer- 
ometric measurements of HD 52961, where the crystalline 
grains seem to be concentrated in the inner disc regions 
(|Deroo et al.ll200^ . The peak position of the 10 /im fea- 
ture in IRAS 06338 points to a lack of amorphous olivine, 
with a preference for amorphous pyroxene and silica, or 
alumino-silicates, in stead. However, both at short and long 
wavelengths we find clear features of forsterite, the crys- 
talline olivine form. HD 52961 shows a 10 /im feature which 
is considerably less broad than other post-AGB disc sources 
(jGielen et al.l l201l[ ) , again pointing to the dominance of 
amorphous pyroxene and silica over amorphous olivine. The 
star also has unusually strong features at 11.3 and 16 //m, 
wavelengths typically associated with forsterite emission. 
The exact shape and strength of these features proved dif- 
ficult to model with current synthetic spectra of forsterite. 

The underabundance of amorphous olivine and weak 
crystalline silicate features in IRAS 06338 and HD 52961 
show that the dust is rather metal-poor. It remains un- 
clear whether this peculiar mineralogy could be linked to 
the extreme chemistry in the stellar photosphere. The pho- 
tosphere for HD 52961 s hows a strong deficien cy of met- 
als, with [Fe/H] = -4.8 (jWaelkens et al.lll991bD . which are 
currently locked in dust grains in the disc. IRAS 06338, 
however, docs not show such extreme depletio n, rather its 
value of [Fe/H] = -1.6 and [Mg/H] = -1.5 (Hr ivnak et all 
I2008D are average valu es found for post-AGB disc sources 
(|De Ruvter et^^m^ . Furthermore, HR 4049 has a simi- 
lar ext reme depletion pattern as HD 52961 (jWaelkens et al.l 
ll991aD . but has a very different infrared spectral signature, 
with features only due to CO2 gas and class B PAHs, and 
no silicate emission. Unfortunately, due to this strong de- 
pletion pattern, it is very difhcult to determine the initial 
metallicity of the central star. 



the class C profile, we would expect more aliphatic bands 
at 6.85 and 7.25 /im, which are (tentatively) observed in 
IRAS 06338 and IRAS 13258. 



7.3. Fullerenes 

Even though the diagnostic value of the Cgo bands in 
IRAS 06338 and HD 52961 is rather hmited, the derived low 
excitation temperatures show that the fullerenes cannot be 
located very close to the central star. Adopting a luminos- 
ity of 5000 L0, a dust temperature of 150 K corresponds 
to a distance of about 240 AU. Unfortunately, neither for 
the PAH or fuUerene features, we can determine their exact 
spatial location. 

Several possible mechanisms for the formation of 
fullerenes i n the C E of evolved stars are discussed by 
ICami et al.l (j2011af ). The cl assic experiments that led to 
the discovery of fullerenes (|Kroto et al.l |1985[ ) show that 
fullerenes self-assemble from collisions between carbon clus- 
ters in a H-poor environment. The presence of H under 
similar conditions inhibits the formation of fullerenes, and 
instead results in the formation of PAHs (Kroto & McK^ 
il988: de Vrie s et al.l 119931 ) . At much higher temperatures 
(> 3 000 K), fullerenes do form in the presence of hydro- 
gen (jJager et al.l l2009h . None of these mechanisms seem 
to correspond to the conditions that could be expected 
in the outflows of evolved stars though. The simulta- 
neo us presence of fullerene s and PAHs in several PNe 
led[G arcia-Hernandez et al.l (|2010[ ) to conclude that PAHs 
and fullerenes form from photoprocessing of hydrogenated 
amorphous carbon (HAG), as was experimentally suggested 
by e.g. Sco tt et al.l ([19971 . However, if this would be the 
fuUerene formation process in PNe, it is hard to under- 
stand why fullerenes are only detected in low-excitation 
PNe, and it would be even more unlikely to invoke such 
a process for the formation of fullerenes in our post- 
AGB sources, where the UV irradiation is considerably 
weaker. One of the more promising avenues is the more re- 
cent study of iMicelotta et al.1 (|2010[ ). that shows fullerenes 
could form from the destruction of PAHs due to shocks. 
IMicelotta et al.l (|2010f ) find that shocks with velocities be- 
tween 75 and 100 km/s are needed to modify PAH molecules 
in the ISM. For several post-AGB disc sources, shocks 
due to strong winds (reaching speeds well above 100 km/s) 
and pulsation (due to the crossing of t he Gepheid insta- 



bility strip) have been measured (e.g . Gillet et all 119891: 
iVa n Winckel"iran [19981 iMaas et al.iT2003t iVan WhM 



20031: iMaas et al.ll2005D . Thus, this scenario could well be 
applicable for the shocked material around these post-AGB 
binaries, where the higher material densities involved could 
be beneficiary to the fuUerene formation. 



7.2. PAHs 

The PAH emission seen in all sample stars can be identi- 
fied as class G (o r BC), dep ending on the underlying sil- 
icate emission (Peetersll201lL and references therein). All 
our sources follow the trend between the effective temper- 
ature of the central star and the central wavelength of the 
7 — 8 fj,m PAH feature (Sloa n et all 120071 [Smolders et"all 
l20ld) . The carriers of the class G fea tures are proposed to 
be aliphatic carbonaceous material ([Sloan et al.ll2007l) . If 
indeed aliphatic carbonaceous material is responsible for 



7.4. Carbonaceous molecules in an O-rich environment 

Given the O-rich classification of IRAS 06338, HD 52961 
and EP Lyr, there must be a process that can free some of 
the locked carbon in the CE to form the observed large car- 
bonaceous molecules. CO molecules can be destroyed due 
to photo-dissociation, but the UV radiation of these sources 
might not be strong enough to dissociate large quantities 
of CO. Alternatively, X-ray radiation, cau sed by shocks, is 
able to break up GO ([Woods et al.ll2005[) . Another possi- 
bility is the Fischer- Tropsch catalysis mechanism, a chem- 
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Table 2. Overvi ew of the different chara cteri stics of post-AGB di sc sources with evidence for carbonaceous molecules from 
the list given bv lDe Ruvter et all (|2006(l and lGielen et alJ (|201l[ ). We list the name, spectral type, effective temperature, 
C/0 ratio, the L/_r/L* ratio, the presence of amorphous and/or crystalline silicates, the detection of PAHs and/or Ceo, 
and the presence of CO2 gas lines in their infrared spectra. 



Name 


Spec. Type 


Teff- (K) 


C/0 


Lir/U 


Am. Sil. 


Cryst. Sil. 


PAHs 


Ceo 


CO2 


EP Lyr 


A4I 


7000 


< 1 


3% 


no / weak 


strong 


class C 


no 


yes 


IRAS 06338 


F3I 


6250 


< 1 


3% 


yes 


weak 


class (B)C 


yes 


yes 


HD 52961 


F6I 


6000 


< 1 


12% 


yes 


weak 


class C 


yes 


yes 


HR4049 


A6I 


7500 


< 1 


25% 


no 


no 


class B 


no 


yes 


IRAS 13258 


F4Ib-G0Ib 


? 


? 


> 50% 


? 


strong 


class C 


no 


no 


Red Rectangle 


FlI 


7500 


< 1 


> 100% 


no 


strong 


class B 


no 


no 



ical reaction process that can form hydrocarbons and H2O 
out of CO and H2. This mechanism has been invoked to 
explain the presence of C-rich mol ecules in other O-rich 
environments an d the solar nebula (jKress fc Tielensll2001l : 
I Jura et ani2006[ ). PAH s have also been found around O-rich 
PNe central stars (e.g. iGuzman-Ramirez et al.ll201ll ). usu- 
ally surrounded by a dense torus. They present a chemical 
model that shows that hydrocarbon chains can form within 
an O-rich environment, but high UV radiation is necessary 
to break up CO, which is not present in our post-AGB disc 
sources. 

Another scenario could be that the central stars have 
become C-rich, but that due to the reaccretion of O-rich 
gas from the circumbinary disc, an opti cally thick layer is 
newly formed on the stellar photosphere (jMathis fc LamersI 
Il992f) . The effects of this accretion process can be seen 
in the photospheric depletion of metals, but it is unclear 
whether the reaccreted gas is rich enough in oxygen, or the 
accretion rate high enough, to explain the apparent O-rich 
stellar photosphere. Furthermore, we ki iow that post-AGB 
stars often undergo strong pulsations ()Van Winckel et aTl 
[1999', 2009), which might easily mix the different layers of 
the stellar photosphere. For IRAS 06338 and HD 52961 pul- 
sation periods of, respectively, 29 and 72 d ays are found 
(jWaelkens et al.lll991bl: [Hrivnak et al.ll2008[ ). 

In the case of IRAS 13258 where the C/0 ratio is not 
known, and the central sources could thus be carbon rich, 
the carbonaceous molecules can be located in a more recent, 
C-rich outflow, while the silicate dust is located in the O- 
rich disc. 



7.5. Relation to other characteristics 

The question remains why C-rich molecules are observed 
in only a few of t he larger sample of ^ 7 p ost-AGB disc 
source s as given in lDe Ruvter et ahl ([200^ and lGielen et al.l 
(|2011h . An overview of the different observational charac- 
teristics of these six sources is given in Table [2J 

So far, most of the post-AGB disc sources where PAH 
emission is seen, also show mid-infrared CO2 gas emission; 
IRAS 13256 and the Red Rectangle being the exception. 
No CO2 gas lines are seen in other sources fro m the larger 
Spitzer sample presented in lGielen et al.l (|201lh . This could 
mean that the CO2 gas is linked to the formation of car- 
bonaceous molecules in an O-rich CE; or that an underlying 
mechanism gives rise to the different observed characteris- 
tics. The presence of CO2 gas could be a by-product of the 
Fischer- Tropsch (FT) catalysis. 



Besides the CO2 and PAH emission, IRAS 06338, 
HD 52961 and EPLyr also all show unusually low L/^/L* 
ratios (see Sect. d]). So far, it is unclear whether the low 
LirI L^, ratio in these sources is a result of a difference in 
disc formation or a signature of disc evolution. The low in- 
frared excess could be due to a lower dust mass in these 
discs, compared to the larger sample. A lower dust-to-gas 
ratio might then make the gas lines more prominent, ex- 
plaining the observed strong CO2 emission in these sources. 
The low infrared excess might point to a disc structure, 
which is more transparent to stellar photons. This would 
then increase the photodissociation rate of CO. 



8. Conclusions 

So far, more than 60 i nfrared spectra o f pos t-AGB 
disc sources, as listed by iDe Ruvter et al.l (|2006t ). have 
been studied, bot h in our ow n Gala xy and the Large 
Magellanic Cloud (IWaters et a l. 1998; Dominik et alll2003l : 
iGielen et al.l l200"8l 2011 ). The spectra are nearly always 
dominated by emission due to O-rich silicate features, with 
a high degree of crystalline dust. Of this large sample, the 
total number of sources where we find evidence for species 
besides silicate dust now comes to six, including the new 
objects discussed in this paper. 

In these six sources the additional features are due 
to O-rich CO2 gas lines and/or carbonaceous molecules, 
most often PAHs. Our detection of Cgo in HD 52961 and 
IRAS 06338 represents the first time these molecules are 
found in binary post-AGB stars. The presence of large aro- 
matic molecules in these environments is very surprising 
and puzzling, given the O-rich nature of the circumbinary 
disc, and possibly the central star itself. 
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Appendix A: Tables 



Table A.l. Best fit parameters deduced from our full spectral fitting. The abundances of small (0.1 /im) and large 
(2.0 nm) grains (S-L) of the various dust species are given as fractions of the total mass, excluding the dust responsible 
for the continuum emission. 





S-L 


Me-FpPvr 
S-L 


Me 

s 


Oliv 
'- L 


S - L 


Silira 
S-L 


Forst 
S-L 


Enst 
S-L 


NaAlSi/iOin 


IRAS 06338 
10 /ttm region A 
10 /ttm region B 
20 /ttm region 


0.0 - 0.0 
0.0 - 0.0 


68.3 - 0.0 
1.2-0.0 


0.0 
0.0 


-0.0 
-0.0 


0.0-0.0 
0.0 - 0.0 


0.0-21.3 
0.0 - 0.0 
0.0 - 16.5 


6.6 - 0.0 
11.6 - 21.1 
15.7-41.3 


1.4 - 2.4 
3.1 - 10.9 
0.0 - 26.5 


51.4 


HD 52961 
10 nm region 
20 /Ltm region 


0.0 - 0.0 


0.0-0.0 


32.0 


- 0.0 


37.6 - 0.0 


0.0 - 10.0 
0.0 - 26.5 


20.4 - 0.0 
17.8 - 40.2 


0.0 - 0.0 
0.0-15.5 




IRAS 13258 
20 /Lim region 












0.0-0.0 


22.8 - 29.6 


9.5 - 38.2 
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